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prove  useful  in  quantifying  the  effects  of  optical  turbulence 
within  the  surface  boundary  layer  on  laser  and  optical 
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ABSTRACT 


Atmospheric  density  fluctuations  induce  phase 
perturbations  that  degrade  the  spatial  coherence  of  a  laser 
beam  propagating  through  the  atmosphere.  These 
degradations  spread  the  laser  beam  and  alter  the  centroid 
and  intensity  profile  stochastically .  Turbulent  conditions 
are  found  at  virtually  all  levels  of  the  atmosphere.  A 
substantial  fraction  of  the  optical  turbulence  along  a 
vertical  path  arises  from  the  heat  flux  between  the 
atmosphere  and  the  Earth’s  surface.  This  type  of 
turbulence  is  typically  within  the  first  100  to  200  meters 
above  the  surface. 

During  this  thesis  research,  a  high  frequency  acoustic 
sounder  was  developed  to  analyze  this  turbulent  layer.  The 
primary  focus  was  the  development  of  the  command  and 
control  software  required  to  coordinate  the  data  collection 
and  reduction.  The  system  was  used  at  two  sites  and  should 
prove  useful  in  quantifying  the  effects  of  optical 
turbulence  within  the  surface  boundary  layer  on  laser  and 
optical  system  performance. 
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A  coherent  laser  beam  propagating  through  the 
atmosphere  is  very  susceptible  to  numerous  turbulence 
inflicted  degradations.  As  electromagnetic  waves  transit 
the  turbulent  regions,  atmospheric  irregularities  randomize 
the  amplitude  and  phase  of  the  wave.  In  order  to  quantify 
the  altitude  dependence  of  the  atmospheric  turbulence,  a 
high  resolution  turbulence  profiler  is  needed.  Acoustic 
echosounders  are  frequently  used  to  detect  and  measure 
atmospheric  density  and  velocity  irregularities  resulting 
from  air  currents,  temperature  inversions,  humidity 
variations,  mechanical  turbulence  and  other  density 
fluctuations . 

Presently  two  atmospheric  optical  parameters,  the 
spatial  coherence  length  ( ro )  and  the  isoplanatic  angle 
(80)  are  measures  of  the  perturbation  of  an  electromagnetic 
wave  propagating  through  the  atmosphere,  and  are  accurately 
measured  by  optical  systems  developed  by  Walters  [Refs.  1 
and  2].  Although  these  systems  measure  a  path  integral  of 
atmospheric  turbulence  with  high  accuracy,  a  major  drawback 
of  each  system  is  that  no  provision  is  available  for 
determining  the  height  of  the  atmospheric  disruptions.  If 
these  disturbances  are  found  to  exist  very  near  the 
surface,  it  may  be  possible  to  negate  their  effect  by 


elevating  the  optical  systems  or  controlling  the  generation 


of  turbulence. 


This  thesis  deals  with  the  design,  construction  and 


implementation  of  a  high  frequency,  acoustic  echosounder 


which  will  accurately  analyze  the  atmospheric  density 


fluctuations  within  approximately  200  meters  of  the 


surface.  As  this  project  is  a  product  of  the  research  and 


efforts  of  two  students,  the  work  was  appropriately 


divided.  My  particular  task  was  to  devise  a  computer 


program  which  would  control  all  aspects  of  the  echosounder 


operation.  Areas  of  particular  interest  were  the 


integration  of  hardware  and  software  within  the  data 


acquistion  system,  developing  plotting  algorithms, 


controlling  the  system  timing,  and  setting  the  input 


parameters  which  determine  the  range  and  sensitivity  of  the 


device.  My  colleage,  LT.  Weingartner  [Ref.  3],  dedicated 


his  efforts  toward  the  actual  design  and  hardware 


development  of  the  echosounder. 


Acoustic  echosounders  have  been  developed  and  in  use 


for  many  years  and  have  proven  to  be  valuable  probes  for 


analyzing  the  structure  and  dynamics  of  the  lower 


atmosphere  [Refs.  4  through  7].  Devices  similar  to  ours 


have  been  used  to  obtain  profiles  of  the  atmospheric 


density  and  temperature  fluctuations  [Refs.  8  and  9],  Our 


device  used  a  high  speed  HP217  computer  to  control  and 


monitor  the  echosounder  which  provided  real  time 


hf’Ti 


m 


information.  Additionally,  we  had  the  ability  to  store  and 


reproduce  the  atmospheric  profile  plots  at  will.  This 


information  was  crucial  in  analysing  the  atmospheric 


turbulence  degradations  on  laser  and  electro-optical 


system’s  performance, 


This  document  shall  address  the  theoretical  background 


of  acoustic  sounder  operation  in  Chapter  II.  Chapter  III 


shall  present  a  synopsis  of  the  system  design  and 


associated  software.  A  summary  of  actual  data  collected  at 


two  sites  is  analyzed  in  Chapter  IV  and  Chapter  V  discusses 


the  conclusions  and  recommendations  found  as  a  result  of 


our  research. 


V.V.tj’, 


Acoustic  echosounders  probe 
transmitting  a  pulse  of  acou 
subsequently  scattered  back  f 


at  a  distance  R  and  frequency  f, 

c  is  the  local  speed  of  sound  (meters/second), 
x  is  the  pulse  length  (seconds), 

A  is  the  aperture  area  of  the  antenna  (meters2), 

R  is  the  range  (meters),  and 

G  is  the  effective-aperture  factor  of  the  antenna. 
Empirically  measuring  or  calculating  the  values  for  all 
other  terms,  one  can  use  this  equation  to  determine 
cro(R.f),  the  scattering  cross  section  per  unit  volume;  that 
is,  the  fraction  of  incident  power  backscattered  per  unit 
distance  into  a  unit  solid  angle  at  a  frequency  f.  Based 
upon  experimental  results,  Tatarski  [Ref.  11]  expresses  the 
acoustic  backscatter  cross  section  per  unit  volume, 
CTO(R.f),  in  the  equation, 

Ct  2 

oo  (  R,  f  )  -  0.00  39  ki  /3  -  , 

To2 

where 

k  -  2k/  A  is  the  incident  acoustic  wavenumber  at 

wavelength  A, 

To  is  the  local  mean  temperature  in  degrees  Kelvin,  and 
Ct 2  is  the  temperature  structure  parameter. 

Combining  this  equation  with  the  echosonde  equation,  one 
obtains  a  volume-averaged  measure  of  Ct 2 . 


11  2  1  Pr 

- - To2  k" 1  /  3 - - - R2exp(2aR) 

0.0039  ErEt  ct  AG  Pr 


Hall  and  Wescott  [Ref.  12]  calculated  a  beam-shape 
compensation  factor  of  0.40  for  a  piston  source  antenna 
with  a  uniformly  illuminated  square  aperture.  This  value 
is  the  same  as  the  effective-aperture  factor,  G,  and  can  be 
substituted  into  the  above  equation.  Approximating  the 
aperture  area  of  the  antenna  to  be  equivalent  to  25  times 
the  aperture  area  of  a  single  speaker  having  a  diameter  of 
7.620  centimeters,  we  get  a  value  of  0.1140  square  meters. 
Combining  these  values  with  the  numerical  constants  in  the 
above  equation,  we  can  simplify  the  equation  for  Ct2. 


1  1  Pr 

Ct2  =  11245  . To2  k-i/3 . —  R2exp(2aR) 

Er  Et  ct  Pt 


The  acousto-electrical  efficiency  factors  for  our 
echosounder  were  measured  with  a  calibrated  microphone  and 
an  anechoic  chamber.  The  calculated  values  for  Er  and  Et 
are  each  0.5.  Substituting  these  values  into  the  above 
equation,  we  can  further  simplify  the  equation  for  Ct2  . 


To2  Pr 

Ct2  =  44980  k-i/3 - - R2  exp(2aR) 


The  above  equation  was  incorporated  into  the  computer 
program  that  provides  the  system  control,  data  acquistion 
and  data  processing  techniques  [Appendix  A].  The  reduced 
data  was  then  used  to  measure  the  temperature  structure 
parameter  as  a  function  of  time  of  day  and  altitude  for 
various  sites. 
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A.  HARDWARE 

Most  of  the  hardware  used  in  this  system  was  standard 
scientific  equipment  and  is  illustrated  in  Figure  1.  The 


Fig.  1.  Echosounder  Equipment  Arrangement 


if* 


pieces  which  were 


special ly 


designed  by  LT  Weingari ner 


I'Ref.  3]  are  the  acoustic  array  and  the  enclosure  for  this 
array.  The  design  considerations  for  the  speakers  and  the 
array  format  are  given  in  Appendix  B  and  the  design  of  the 
enclosure  to  house  this  array  is  outlined  in  Appendix  C. 
Aside  from  these  two  pieces  of  equipment,  there  are 
primarily  six  other  components  which  comprise  this 
echosounder.  Below  is  a  brief  description  of  each  of  these 
additional  components. 

1 .  HP200  Series  Computer 

The  Hewlett-Packard  (HP)  200  Series  Computer 
includes  a  20  megabyte  hard  disk  drive,  a  floppy  disk  drive 
and  an  associated  printer  and  monitor.  The  HP200  Series 
Computer  was  the  central  control  component  for  the  entire 
echosounder  arrangement.  The  computer  used  was  an  HP21/ 
programmed  in  Basic  3.0  and  equipped  with  an  Infotek  80203 
Basic  Compi Ler  and  an  Infotek  FP210  Floating  Point 
Accelerator  to  enhance  the  speed  of  program  execution.  The 
program  defined  all  parameters  for  the  HP3314A  Function 
Generator  as  well  as  execute  the  trigger  command  wh:  ,;h 
produced  the  echosounder  transmitted  acoustic  puLse.  The 
computer  also  received  the  data  from  the  acoustic  array  v i a 
the  pre-amplifier,  bandpass  filter  and  analog  to  digital 
river t--r.  The  computer  then  conducted  the  data  reduction 
rout- ir.es  to  produce  and  display  a  high  resolution 
a tmosphe r ic  profile. 


lb 


2 .  HP3314A  Function  Generator 

The  Hewlett-Packard  (HP)  3314A  Function  Generator 
was  a  multimode  function  generator  capable  of  providing 
sine,  square  and  triangular  wave  functions  as  well  as  any 
desired  waveform  ranging  in  frequency  from  0.001  Hertz  to 
19.999  Megahertz.  The  HP3314A  Function  Generator  was  used 
to  supply  a  pulse  of  an  integer  number  (usually  100)  of 
sinusoidal  cycles  of  constant  amplitude  to  the  QSC  Model 
1700  Audio  Amplifier.  A  constant  frequency  setting  of  5000 
Hertz  was  used  for  all  data  runs. 

3 .  QSC  Model  1700  Audio  Amplifier 

The  QSC  Model  1700  Audio  Amplifier  is  a  high  power 
amplifier  which  can  supply  350  watts  over  an  8  ohm  load. 
This  amplifier  was  used  to  boost  the  output  voltage  of  the 
function  generator  by  a  factor  of  20  from  1.5  volts  to  30 
volts.  The  power  supplied  to  each  speaker  in  the  array 
during  operation  was  then  37.5  watls. 

4  -  Pre- amplifier 

The  pre-amplifier  was  designed  and  constructed  by 
Walters  [Ref.  13]  to  supply  a  gain  of  1000  to  the  returned 
signal  and  isolate  the  data  acquistion  components  from  the 
transmitted  pulse.  The  pre -ampl i f ier  may  be  thought  of  as 
a  safety  and  switching  mechanism  for  the  system.  An  LT 
1037  and  an  LT  1007  Operational  Amplifier  were  selected  for 


use  in  the  pre-amplifier  based  upon  their  low  noise 
properties  which  were  evaluated  by  physically  incorporating 


them  into  the  device  and  measuring  their  noise  and  gain 


characteristics  with  an  HP3561  Spectrum  Analyzer.  The 
diodes  used  are  IN  4000  Series  Rectifiers  capable  of 
carrying  one  ampere.  These  rectifiers  serve  to  isolate  the 
power  amplifier  from  the  amplifier  input  stage  and  to  limit 
the  voltage  applied  to  the  operational  amplifiers.  The 
electrical  diagram  of  this  device  is  outlined  Figure  2. 


Fig.  2.  Schematic  of  Pre-amplifier 
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The  Rockland  Wavetek  Model  852  Filter  operated  as  a 
48  db  per  octave  bandpass  filter  to  suppress  the  broadband 
noise  of  the  system.  High  and  low  bandpass  settings  of 
5500  Hertz  and  4500  Hertz  respectively  were  used.  The 
filter  response  at  these  settings  is  illustrated  in  Figure 
3. 


START  1  000  Hz  RANGES  -5  ctBV  STATUS:  PAUSED 

A:MAG  RMS: 50 


START:  l  080  Hz  8U:  477.43  Hz  STOP:  31  000  Hz 

X:  3000  Hz  Ys-28.24  dBV  _ 


Fig.  3.  Bandpass  Filter  Response 

The  three  db  (half  power)  bandwidth  is  approximately  1  KHz 
wide.  Although  it  would  be  desirable  to  reduce  the 


bandpass  to  between  50  and  100  Hertz  without  impairing 
signal  response,  the  Rockland  Wavetek  Model  852  Filter  was 


the  best  filter  available.  It  performed  adequately  in 
spite  of  the  large  lKHz  bandwidth. 

6 .  Infotek  AD2QQ  Analog  to  Digital  Converter 

An  Infotek  AD200  12  bit  Analog  to  Digital  Converter 
was  used  to  digitize  the  signal  voltage  for  the  computer 
system  at  a  12.5  to  20  KHz  sample  rate. 

B .  SOFTWARE 

The  software  and  HP200  Series  Computer  system  were 
responsible  for  controlling  and  monitoring  every 
operational  phase  of  the  hardware  components.  The  program 
which  accomplished  this  task  was  entitled  "ACRDR"  and  is 
listed  in  Appendix  A.  This  program  was  modeled  after  a 
program  written  by  Walters  [Ref.  13]  but  each  program 
performs  a  distinctly  different  computational  task.  This 
program  was  written  in  HP  Basic  3.0  and  was  compiled  by  an 
Infotek  BC203  Basic  Compiler  to  enhance  the  speed  of 
execution . 

The  program  "ACRDR"  is  easily  broken  down  into  a  number 
of  blocks  and  subroutines  which  performed  specific 
operations.  These  sections  are  outlined  in  a  flowchart 
(Fig.  4).  Each  block  is  straightforward  in  its  purpose, 
and  the  program  is  designed  to  be  as  helpful  to  the  user  as 
possible.  As  a  prologue  to  the  actual  program  code,  there 
is  a  listing  of  all  the  program  variables  with  a  short 
description  of  their  use.  Such  a  listing  familiarizes  a 
user  with  the  computations  to  be  made  and  also  provides 


19 


quick  references  for  any  future  modifications  to  be 
implemented . 


Fig.  4.  Flowchart  of  Fundamental  Program  Operations 

The  initialization  procedures  which  set  all  the 
parameters  necessary  for  data  collection  follow  the 
variable  definition  section.  These  parameters  control  such 
things  as  the  contrast  between  background  and  return  signal 
and  the  setting  of  the  computer’s  internal  clock.  All 


options  of  the  HP3314A  Function  Generator  are  decided  by 
the  user  at  this  phase.  Another  initialization  involves 
preparing  the  Infotek  AD200  Analog  to  Digital  Converter  for 
operation.  Finally,  internal  arrays  are  dimensioned, 
integer  variables  are  defined  and  the  computer  function 
keys  are  redefined  to  suit  the  "ACRDR"  program. 

Now  that  the  system  is  ready  to  operate,  certain 
quantities  must  be  determined  based  upon  atmospheric 
conditions  at  the  time  of  data  collection.  From  the 
ambient  temperature,  a  speed  of  sound  in  air  is  determined 
by  the  equation  in  Kinsler,  et  al .  [Ref.  14:  p.  106] 

C  -  Co  [  1  +  (T/273)  }i  /2  , 

where 

C  is  the  local  speed  of  sound, 

Co  is  the  speed  of  sound  in  air  at  0"  Celsius  (331  6 
m/s ) .  and 

T  is  the  local  temperature  in  degrees  Celsius. 

Using  this  speed  with  a  user  decided  sampling  frequency  and 
number  of  data  points,  the  maximum  range  of  the  traces  is 
determined  by  the  relation, 

RANGE  =  C/2  •  1/ ( SAMPLE  FREQUENCY)  •  N  , 

whe  re 

RANGE  is  the  range  of  the  echosounder  in  meters, 

SAMPLE  FREQUENCY  is  the  user  desired  rate  of  data 
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acquisition,  and 

N  is  the  number  of  data  points  per  trace. 

Typically,  the  values  of  many  of  the  parameters  are  not 
varied.  The  value  of  N  above  is  almost  always  taken  to  be 
the  allowed  maximum  of  16301  data  points  and  the  speed  of 
sound  can  be  “afely  estimated  to  be  340  m/s  at  room 
temperature.  By  using  a  typical  sample  frequency  of  20000 
Hz,  a  range  of  approximately  135  meters  is  obtained. 

Another  parameter  determined  is  the  atmospheric 
attenuation  coefficient.  This  calculation  is  done  in  a 
subroutine  obtained  from  Reference  15  and  is  necessary  in 


this  program  for  the  calculation  of  Ct 2 . 


Other 


computations  made  prior  to  program  execution  include  the 
pulse  length,  wavenumber,  and  the  D-C  offset  of  the 
equipment . 

At  this  point,  the  computer  is  finally  ready  for  data 
collection.  The  screen  setup  displays  a  distance  versus 
time  plot  and  the  internal  clock  of  the  computer  is 
synchronized  with  the  pulse  from  the  HP3314A  Function 
Generator.  Following  transmission  of  the  pulse,  the  return 
signal  is  received,  digitized,  stored  in  an  array  and  data 
reduction  commences.  A  block  averaging  technique  is  used 
in  which  a  block  of  data  points  the  size  of  the  number  of 
emitted  cycles  is  summed  and  averaged.  Returns  are  plotted 
as  darkened  areas  on  the  echosounder  traces  with  the 
intensity  of  the  darkened  area  linearly  related  to  the 
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magnitude  of  the  return  signal.  This  same  procedure  is 
used  in  the  computation  of  Or2  with  some  very  important 
differences.  Aside  from  simply  correcting  for  7rra 
spherical  divergence.  Ct 2  is  corrected  for  electronic  gain, 
the  ratio  of  power  returned  to  power  transmitted,  the 
efficiency  of  the  speakers  as  transmitters  and  receivers, 
the  area  of  the  speaker  array,  the  atmospheric  attenuation, 
the  pulse  length,  the  temperature,  the  scattering  cross 
section  per  unit  volume  at  a  specific  range  and  frequency 
and  finally  the  effective  aperture  factor  of  the  antenna. 
Additionally ,  Ct2 (Range)  is  averaged  for  a  particular 
altitude  over  15  minute  intervals.  Finally,  after  each 
pulse  is  reduced,  a  corresponding  mean  square  noise  level 
is  determined.  This  is  done  by  averaging  each  block 
average  at  maximum  range  until  at  least  ten  values  have 
been  used  in  the  average.  An  upper  limit  of  five  over  the 
average  noise  figure  (this  corresponds  to  voltage 
fluctuations  on  the  order  of  ID-7  volts)  is  set  on  the 
routine  to  avoid  averaging  any  strong  return  signals  or 
anomal ies  such  as  passing  aircraft.  After  the  ten  values 
are  averaged,  every  subsequent  pulse  is  averaged  into  all 
the  preceding  noise  levels  and  removed  from  each  subsequent 
return  signal. 

Ultimately,  after  each  15  minute  interval,  or  at  the 
users  request,  the  display  terminal  image  is  printed.  At 
this  time  the  Ct 2  computations  are  conducted  and  plotted  on 
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California.  This  site  is  located  atop  Mt .  Hamilton  at  an 
altitude  of  approximately  5700  feet  and  nearly  20  miles 
inland  from  the  coast. 

These  two  data  collection  sites  represent  areas  of 
differing  atmospheric  air  pressures,  water  vapor  pressures, 
local  temperature  ranges,  and  local  wind  velocity  ranges. 
These  characteristics  all  play  important  roles  in  effecting 
the  local  atmospheric  turbulent  conditions  and  thereby  the 
atmospheric  structure  parameter,  Ct2. 

In  addition  to  collecting  echosounder  data  at  Lick 
Observatory,  simultaneous  measurements  of  the  isoplanatic 
angle  (So)  and  spatial  coherence  length  ( ro )  were  made  with 
systems  developed  by  Walters  [Refs.  1  and  2],  A  basic 
knowledge  of  these  two  systems  is  necessary  to  understand 
the  correlation  procedures  made.  The  isoplanatic  angle 
(0o)  is  primarily  an  upper  atmospheric  measurement  which 
indicates  atmospheric  disruptions  at  a  range  of  2  to  15 
kilometers.  The  spatial  coherence  length  ( ro )  is  a  measure 
of  the  effects  of  the  entire  atmospheric  blanket  on 
coherent  light  transmission.  A  close  comparison  of  all 
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the  low  altitude  boundary  layer  contribution  is 
sufficiently  small.  This  strong  correlation  helps  to 
reinforce  the  overall  description  of  the  atmosphere  at  the 
time  of  data  collection.  Figures  5  and  6  graphically 
illustrate  the  atmospheric  measurements  made  at  Mt.  Wilson 
on  2  April  1987  by  Walters. 
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Fig.  6.  Spatial  Coherence  Length  Measurements,  Mt.  Wilson 


The  strong  correlation  between  the  isoplanatic  angle 
(8o)  and  the  spatial  coherence  length  (ro )  is  especially 
evident  between  the  hours  of  0700  and  1300  universal  time. 
The  close  tracking  of  these  two  measurements  during  this 
time  interval  indicate  that  the  upper  atmospheric 
conditions,  as  measured  by  0o ,  are  dominating  the  entire 
atmospheric  profile  as  measured  by  ro .  Unfortunately,  the 
existence  or  non-existence  of  any  turbulent  surface  effects 
cannot  be  ascertained  by  the  employment  of  the  above  two 


systems  alone.  However,  employment  of  these  two  systems 


together  with  the  acoustic  echosounder  should  enable  us  to 
produce  a  complete  atmospheric  profile  with  strong 
correlation  between  all  three  atmospheric  measurements. 

On  9  and  10  April  1987,  all  three  systems  were  operated 
at  Lick  Observatory.  Again,  a  good  correlation  between  the 
isoplanatic  angle  (9o )  and  the  spatial  coherence  length 
( ro )  measurements  was  noted.  In  addition,  a  strong 
correlation  between  the  spatial  coherence  length  ( ro )  and 
echosounder  measurements  was  present.  A  comparison  of  the 
atmospheric  data  in  Figures  7  and  8  shows  a  good 
correlation  between  the  two  parameters  especially  during 
the  0930  to  1130  time  interval  on  9  April  1987.  However, 
during  subsequent  hours  the  isoplanatic  angle  ( 9o ) 
measurements  remain  high  (=*12  i_irad)  indicating  relatively 
calm  turbulent  conditions  in  the  upper  atmosphere  while  the 
spatial  coherence  length  ( ro )  values  drop  sharply  after 
12:00  Universal  Time  indicating  dominant  and  increasing 
lower  atmospheric  turbulence.  This  increase  in  the  low 
level  turbulence  should  be  evident  in  the  echosounder  data 
commencing  around  1200  universal  time  (0400  local  standard 
time)  on  9  April  1987.  A  comparison  of  the  echosounder 
data  in  Figures  9  through  11  illustrates  this  increase  in 
the  local  surface  turbulence. 


Fig.  7.  Isoplanatic  Angle  Measurements,  Lick  Observatory 


Fig.  8.  Spatial  Coherence  Length  Measurements 
Lick  Observatory 
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Fig.  9.  Echosounder  Trace,  Lick  Observatory 
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Fig.  11.  Echosounder  Trace,  Lick  Observatory 

The  strong  correlation  between  the  echosounder  data  around 
4:00  to  5:00  Standard  Time  and  the  ro  measurements  around 
12:00  Universal  Time  combined  with  the  lack  of  correlation 
between  the  0o  and  ro  measurements  indicate  that  the  lower 
atmospheric  and  surface  turbulence  are  dominating  the 
atmospheric  profile  during  this  time  period. 

Data  collected  on  10  April  1987,  again  illustrate  the 
strong  correlation  between  the  three  atmospheric 
measurements  made.  A  comparison  of  the  9o  and  ro 


measurements  during  the  time  interval  of  0700  and  1300 
Universal  Time  indicate  steady  turbulent  conditions  in  the 
upper  atmosphere  and  greatly  varying  turbulent  conditions 


at  lower  atmospheric  levels.  This  is  evident  in  Figures  12 
and  13  by  Lhe  consistent  values  of  00  during  the  time 
period  compared  with  the  steady  increase  and  eventual 
decline  of  the  ro  values  during  the  same  time  interval. 
The  trace  variations  in  Figure  13  during  the  hours  of  0800 
and  1200  are  indicative  of  a  period  of  decreasing  lower 
atmospheric  or  surface  turbulence  followed  by  the  onset  of 
an  increasingly  turbulent  period  around  1200.  This 
turbulent  trend  is  strongly  supported  by  the  echosounder 
data  in  Figures  14  through  21. 


Fig.  12.  Isoplanatic  Angle  Measurements,  Lick  Observatory 
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Fig.  21.  Echosounder  Trace,  Lick  Observatory 

The  large  200mm  coherence  lengths  around  11:00  Universal 
Time  are  consistent  with  the  low  turbulence  evident  in  the 
echosounder  profiles  around  3:00  Standard  Time  followed  by 
a  pre-dawn  increase  in  the  surface  turbulence. 

Only  acoustic  sounder  atmospheric  measurements  were 
made  on  the  upper  roof  of  Spanagel  Hall  at  the  Naval 
Postgraduate  School.  Data  runs  on  26  and  27  April  1987, 
included  both  echosounder  data  and  the  associated  Cl 2 
plots.  A  representative  sample  of  the  data  collected 
during  this  period  is  included  in  Appendix  D.  Of 


particular  interest  is  the  data  in  Figure  22,  which  shows 
the  maritime  boundary  inversion  layer  at  about  100m  and 
being  perturbed  by  convective  plumes  at  lower  altitudes. 


A  high  frequency  (5KHz)  phased  array  echosounder 


constructed  to  measure  low  level  turbulence  appears  to  work 


within  the  10  to  135  meter  altitude  range.  When  using  the 


device,  detailed  profiles  of  the  short  range  atmospheric 


density  fluctuations  were  obtained.  The  profiles  were 


found  to  correlate  very  well  with  the  measurements  of  the 


isoplanatic  angle  (80)  and  the  spatial  coherence  length 


( ro )  during  periods  of  simultaneous  operation.  This  short 


range  echosounder,  when  used  in  conjunction  with  the  other 


atmospheric  measuring  devices,  is  an  invaluable  tool.  The 


graphical  output  provides  us  with  a  more  complete 


description  of  the  atmosphere  that  can  be  used  to  calculate 


Cn 2  ,  the  atmospheric  index  of  refraction  structure 


parameter . 


Further  software  development  such  as  the  incorporation 


of  Fast  Fourier  Transform  routines  into  the  echosounder 


program  will  determine  the  radial  velocity  profile  of  the 


return  signal.  Such  an  improvement  will  determine  the 


velocity  of  the  probed  air  masses  and  plot  these 


echosounder  traces  as  a  function  of  color  intensity. 


Other  areas  of  further  reasearch  include  the  design  and 


testing  of  different  array  patterns.  One  such  pattern,  the 


hexagonal  array,  is  already  undergoing  tests.  As  the 


speaker  array  evolves,  it  is  also  evident  that  the 
enclosure  must  follow  suit  to  accomodate  the  new  array 
design  and  associated  beam  patterns.  Further  improvements 
in  software  routines  are  also  inevitable.  It  is  always 
desirable  to  store  on  a  disc  all  data  collected  at  a  site. 
Presently,  the  floppy  disc  capacity  is  inadequate  for  the 
storage  of  more  than  two  hours  of  data.  Furthermore,  if 
FFT  and  Doppler  routines  are  to  be  added,  the  computational 
speed  of  the  computer  will  inhibit  the  pulse  repetition 
rate.  A  special  data  acquisition  technique  called  "Direct 
Memory  Access"  may  then  have  to  be  introduced  into  the  HP 
computer  to  allow  simultaneous  data  collection  and 
processing . 

Finally,  this  product  may  be  used  in  various 
applications  not  already  addressed.  One  such  use  may  be  to 
measure  windshear  at  airports.  Other  researchers  have 
expressed  interest  in  using  the  echosounder  to  measure 
arctic  atmospheric  conditions  during  meteorological 
surveys.  Assuredly,  as  the  device  evolves  into  a  more 
compact  and  highly  efficient  instrument,  its  range  of 
application  will  continue  to  expand. 
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The  computer  program  receives  information  from  an  acoustic 
array  through  an  A-D  converter.  This  information  is  the 
returned  signal  of  an  acoustic  pulse  as  it  passes  through 
the  atmosphere.  The  data  is  then  used  to  display  the  return 
intensity  with  distance  as  a  function  of  time.  This  is  a 
short  range  device  (from  150m  to  280m). 


LIST  OF  VARIABLES 


counters  for  loops 


Rec_num  4  Nrec  -  real  and  integer  representation  of  the 

record  number  for  storage 


Plotnum  -  counter  used  to  insert  form  feeds  between  plots 


ui sc _ add'  e335  -  storage  location  of  data  fii; 


1 1  eS 


lie  name  u5Su  to  store  data 


string  and  real  reor ssentat  ion  of  the  number  of 


Cyc.as  and  the  number  of  point3  used  in  the 
computation  of  the  block,  average 


Point(»)  -  real  number  representing  the  noise  and  range 
corrected  average  over  given  number  of  points 


Dat (  *  ) 


array  of  a-d  converter  output  after  sampling 
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(SF 


420 

!  02<  * ) 

array  of  reduced  and  averaged  data  including 

430 

! 

offset,  correction  for  noise,  and  range 

440 

! 

correct  ions 

450 

! 

460 

!  Hrs  - 

the  integer  hours 

470 

! 

403 

!  Min  - 

the  number  of  minutes 

490 

! 

500 

!  Qtrhr 

& 

Qtrmin  -  keeps  track  of  the  passage  of  each 

513 

! 

15  minute  interval 

523 

! 

533 

!  Numl 

- 

counter  for  computation  of  block  averages 

540 

! 

550 

!  Num2 

- 

counter  for  the  number  of  block  averages  made 

5S3 

| 

570 

!  M  - 

the  counter  used  in  the  plot  label  routines 

503 

I 

530 

!  Timcnt 

time  interval  between  data  samples  in 

600 

I 

nanoseconds 

610 

! 

623 

! 

633 

!  Kmax 

- 

total  number  of  block  averages  computed 

640 

! 

650 

!  Freq$ 

& 

Freq  -  the  frequency  input  of  the  HP3314A  in  kiloher 

663 

! 

670 

!  Zone® 

- 

the  appropriate  time  zone  the  operator  desires 

603 

! 

690 

!  Sites 

- 

the  name  of  the  appropriate  site  of  data  collection 

700 

t 

713 

1  2ave_plt 

&  Cntrl  -  on/off  toggles  used  to  determine  whether 

723 

! 

a  particular  run  will  be  saved  to  a  disc 

I  w>  44 

I 

743 

!  Nplot 

- 

the  number  of  increment  along  the  entire 

753 

! 

horizontal  axis 

760 

! 

773 

!  Offset 

the  computed  0  -C  offset  for  the  system  prior 

"iC.n 
(  U  <J 

1 

to  data  acquisition 

733 

I 

0C3 

*  No  i  3  s 

- 

a  running  total  of  the  noise  accumulated  at  the 

n  * 

U  i  u 

! 

rta  A  1  f'lUri  1  31  iy  C  i  u  j  u  w  1  4,  :  :  U  <2  /  .  w  C  £  I'l  G  v  C  i.  he 

C23 

! 

d  v  ci  dye  l  \ 0  i  5  c 

n-rn 

O  U'O 

n  *  r* 

U  -♦'C' 

!  Limit 

- 

the  upper  bound  on  the  noise  figure  used  to  insure 

853 

that  large  returns  are  not  included  in  the 

863 

noise  computations 

370 

800 

!  K 1  - 

a 

counter  of  the  number  of  traces  used  to  determine 

890 

I 

the  noise  average 

310 

1 

Samave 

a  sample  average  used  in  the  computation  of  the 

320 

! 

D-C  offset 

320 

1 

340 

Of f_ave 

a  toggle  used  to  inhibit  the  firing  of  the 

950 

HP3314A  during  0-C  offset  computations 

9S0 

370 

Count® 

u  Cnt  -  the  string  and  real  representation  of 

980 

1 

the  number  of  data  points  to  be  taken 

990 

l 

1000 

! 

Samfreq 

the  user  input  sampling  frequency  desired 

1010 

i 

1020 

! 

Tninc 

1  over  the  sampling  frequency!  the  time  between 

1030 

1 

samples 

1040 

! 

1050 

1 

T  ine$ 

the  string  required  by  the  A-D  converter  to 

1050 

I 

sample  at  the  desired  rate 

1070 

! 

1030 

i 

COkelvin  -  the  speed  of  sound  at  0  degrees  celcius  in  m/s 

1090 

! 

1  100 

! 

Temp 

temperature  of  the  surroundings  in  degrees  celcius 

1110 

! 

1120 

1 

Spd 

the  speed  of  sound  in  air  computed  for  the  input 

1130 

| 

temperature 

1  140 

! 

1  150 

! 

Rspd  - 

the  relative  speed  of  sound  that  the  echosounder 

1180 

! 

sees  which  is  half  the  computed  speed 

1  170 

! 

1180 

! 

Rdist 

the  distance  traveled  in  one  time  increment 

1190 

! 

1200 

! 

Maxrng 

the  maximum  range  of  the  echosounder  found 

1210 

I 

by  multiplying  the  distance  per  time  increment 

t  nn<% 

1  4.4.  V 

1 

1 

by  the  number  of  time  increments 

1  twu 

1240 

I 

Exrng 

the  range  rounded  to  the  nearest  value  evenly 

1250 

! 

divided  by  15;  used  solely  for  plotting 

1280 

1  27(3 

I 

I 

purposes 

*  nr*r% 

I  4.WU 

I 

Far 

used  to  estimate  the  far  field;  this  value  keeps  the 

1220 

! 

range  correction  from  being  applied  to  data  very  near 

1  r\  r% 

l  .uu 

«  ?  •  s 
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to  the  source 

t  w  '  \J 

1  T  •*»  r* 

I  1  v  1 

the  user  input  intensity  level  divisor;  this  value 

1  771 

1 

| 

sets  the  screen  contrast  in  data  return 

1  ■?  *  f* 

1  j'tt; 

| 

1  350 

! 

3c a le _ y 

the  ratio  of  exrng  to  maxrng;  this  value  is  used 

1380 

1  370 

! 

| 

to  keep  the  plotted  range  of  data  accurate 

1330 

i 

Npoint  u 

Npointl  &  Npoint2  -  variables  used  to  keep  track  o 

1390 

l 

time  passage  between  pulses 

1  400 

! 

variables  used  to  keep  track  of 
tine  passage  between  pulses 
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1423 
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T2  &  T1  &  T0  -  used  to  synchronize  data  collection  with  the 
clock 


!  TIMEDATE 


the  internal  clock  of  the  computer 


Ovd  -  divisor  of  the  data  average;  used  because  computer 
multiplication  is  faster  than  division 

R  -  the  range  of  a  block  of  data  samples 

Run_ave  -  the  running  sum  of  the  block  samples 

X  -  the  value  of  the  data  point  less  the  D-C  offset 

Ns  -  the  final  average  of  the  block'  of  data  points;  also 
used  in  the  computation  of  the  noise  figure 

Corr  -  the  noise  correction  applied  to  the  data  samples 

Time  -  the  horizontal  position  of  the  trace  on  the  plot 

Timedist  -  the  horizontal  width  of  the  trace  on  the  plot 

£nd_pt  -  the  last  point  of  the  plot  vertically  taking 
into  account  the  scaling  factor  scale_y 

Inc  -  The  vertical  increment  along  the  plot 

Z  -  The  final  reduced  data  points  which  are  output  on  the  plot 

Dis  -  The  vertical  height  of  the  trace  on  the  plot 

Ap$ ,Frq$ ,Nm$ ,En$,VoS ,Hz$  -  strings  needed  to  set  the  HP3314A 

Amp$  -  userv  input  amplitude  for  the  HP3314A 

Pls_lng  -  the  pulse  length  of  the  burst 

C t ( »  5  -  the  atmospheric  temperature  structure  parameter 

Y  t  ft  c  f.  *v'  u  c  *  0  *  i"  5C  5  3  <_*  —  3  v_i  ill  gwripu  •  *  ny  w  w  »  *  . 


i  n  t  r\ 

\ 

1  v/  —  \ 

1  fpi 

the 

vSrt  iCdl  p03iti0“i 

on  ths  Ct  plot 

t  n  4  f% 

1 

i 

1053 

1  Pnt 

the 

horizontal  position  on  the  Ct  plot 

1363 

l 

J 

\  872 

!  K3 

the 

wave  number  to  the 

1 /3rd  power 

1  nnr* 

1  oow 

1 

1  G  GO 

!  A  - 

The  area  of  the  receiver 

( array  ) 

1320 

l 

45 


1313 

! 

G  - 

The  effective  aperature  factor 

1923 

! 

• 

1933 

t 

Er  - 

efficiency  of  coversion  of  acoustical  power 

1943 

! 

to  electrical  power  on  the  recieve  side 

1353 

! 

1363 

! 

Et  - 

the  efficiency  of  conversion  of  electrical 

1373 

I 

power  to  acoustical  power  on  the  transmit  side 

1  3S3 

t 

1333 

! 

I 

■*-> 

CL 

The  computed  transmitted  power  to  the 

*7  r%  r*f\ 

U  1/1111 

! 

acoustic  array 

n/%  i  r% 
uu  i  u 

! 

-i  run  r% 
LULti 

! 

Gain 

the  electronic  gain  of  the  equipment 

2033 

1 

2340 

! 

Z  imp 

the  speaker  array  impedence 

2053 

I 

2363 

! 

2073 

| 

2333  OPTION  BASE  1 


initialize  the  arrays  &  set  dimensions 
declare  all  integer  variables 


DIN  Disc_addressSC 20 ] ,FileSC33  3 .Point ( 330 > ,Ct ( 203  > 
INTEGER  I .Hr .Plotnum ,Print_key  ,Num  ,Numl  ,Num2  ,N  ,K 
INTEGER  Rec_num ,Kmax 
INTEGER  02(300  ),0at<lG3«n  BUFFER 


initialization  rout  ines ....  set  time,  set  HP3314A  Function 
Generator 


2203  Rstrt:CALL  Freq_init ' Freq* ,N$ ) 
CALL  Init_ad230 
CALL  Set_t ime< Zone*  ) 

INPUT  "SITE  NANE  “.Sites 

v  1 1  *n 


2263 

C  U  I  u 
n  n  33 

»  J  a 

^  XJ 


2330  ! 

2133  I 
25  13  ! 

2123  ! 

2133  ! 

2143  ! 

2153 
2160 
2173 
2180 
2130  I 
2233  ! 

ni  i  m  I 
C.L.  \  11  I 

2223  I 

I 

2243  ! 


2343 

2303 

2363 


2333 

2333 


!  keyboard  3et_up  -  set5  labels  on  the  C  G  rap  Li  t  Si*  fufiCt  i  G  H 
1  keyS 

I 
I 

OUTPUT  KBD; “SCRATCH  KEYE“ ;  !  CLEAR  THE  KEYS 

CONTROL  2  ,2; 1 
ON  KEY  1  LA8EL  “NEW 


TENP “  GOTO  Speed 


ON  KEY  2  LABEL  * INTENS .  FACTOR*  GOTO  Lvl 
ON  KEY  4  LABEL  ‘PRINT  TRACE"  GOTO  Prt_dmp 

ON  KEY  5  LABEL  ‘SAVE  NXT  PLOT*  GOTO  Svpt 
ON  KEY  S  LABEL  “CHANGE  SMPL  FRO"  GOTO  Sfrq 
ON  KEY  7  LABEL  ‘RESTART*  GOTO  Rstrt 
ON  KEY  3  LABEL  ‘QUIT*  GOTO  Quit 


Set  constants 


0i3C_addres3$»“ : ,730 ,1  ,3" 

Maxrec-723 
Np  lot*900 
Plotnum* l 
Of  f set»0 
GINIT 
Noise-O 
Limit*! 3 
Y-3 
Kl-0 
Samave-0 
A-. I  143 
G-.4 
Er*.49S 
Et-.43G 
Gain*731 2S"2 
Zimp*»l  2.12 
! 

! 

!  Input  and  calculation  of  terns  used  by  the  A-D  converter  and 
!  the  computation  of  range 


PRINT  “Enter  the  number  of  data  points  desired  (max  15301)“ 
PRINT  ‘Use  increments  of  the  number  of  cycles  then  add  1“ 
PRINT  “ ex :  (150  cycles  X  133  data  points /cycle) ^1*15331* 
PRINT  "RECOMMEND  USING  MAX  UAL'JE  OF  15301“ 

INPUT  Count! 

I 

'  uutpiit  pwwfiT  CwnuutSd  a55j-'ici 

I  A  -  _  1  .  £  .  -  -  _  x.  A  .  _  —  _ 

nonpu  i  *  Ci  :Cv 

i 

Pt-(EtM30"2  ) )/(  2»2imp  ) 


23S3  Sfrq:PRINT  "Input  the  sampling  frequency  desired.  This  frequency  will" 
2353  ‘  -PRINT  "determine  the  range  of  the  sounder.  Examples  are'." 


2373  PRINT  "  SAM.  FREQ.  -  12500 . RANGE  -  225  M" 

2330  PRINT  “  SAM.  FREQ.  -  23033 . RANGE  -  135  M" 


2330  INPUT  "ENTER  THE  SAMPLE  FREQUENCY"  .Samfreq 


2933  GCLEAR 

2910  Tminc- 1 . /Samfreq 

2923  Timcnt- !  000*  INK  Tminc/ 1  .E-S  ) 

2923  Time*- "TIME  “&UAL$< Timcnt  ) 

2943  Cnt-UAL<  Counts ) 

2953  C3kelvin-331 .5 

2953  Num-UAL(NS) 

2973  Freq-UAL(FreqS )» 1 300 

2933  I 
2990  ! 

3333  !  Computation  of  the  Q-C  offset  prior  to  program  run 

~rr%  ♦  r%  | 

<j  i  %}  ! 


3023 

3343 

3053 


Jt/  O  <J 

3373 


3393 
3133 
31  13 
3123 
3133 
3143 
3153 
3153 
3173 
3183 
3190 
3233 


Off.ave-l 

Offset-3 

FOR  1-1  TO  13 

CALL  Read_ad230( Oat ( » ! .Counts .TimeS ,0f f_ave ) 
FOR  K-1  ’0  ( Cnt-1  ) 

Samave-Samave+Dat  ( K ) 

NEXT  l< 

PRINT  “COMPUTING  D-C  0FF5ET" 

NEXT  I 
Of  f_ave-0 

Of f set-Samave/( 10*(Cnt-t  ) ) 

PRINT  “OFFSET  IS  :“  .Offset 


i 


I  computation  of  the  speed  of  sound  at  a  given  temp  and  the 
!  range  of  detection  of  the  device 


I 


U  L. 


U  -TV 

3243 

3253 


3213  Speed! INPUT  “Enter  the  temperature  (Celsius)  *  .Temp 
Spd-C3keIvin»(SQR( 1 +< Temp/273  ) ) ) 

Rspd-Spd/2 
Rdist-Tminc»Rspd 
Maxrng-Cnt  »Rdist 

Exrng»( ( INT( Maxrng/33 ) )+ . 5  )  »33 . 3 
Far-I ./( Tminc *Spd  ) 

Lambda-Spd/Freq 
K 3  —  < ( 2  »P I ) /Lambda  ) A ( 1 / j  / 

P 1  s  _  1  r.  g  —  (  R  s  P  d  *  N  u  m  )  /  F  r  e  q 
CALL  At t er.ua t  ior,(  Freq  .Temp  , Attar,  i 


jLtfd 
“7  n  r% 


:g3 


3343 

3353 

3353 

3373 

3330 

3333 

3433 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


“Enter  the  relative  intensity  division  level.  Thi, 
“value  is  used  to  determine  the  plot  intensity  by 
“dividing  the  block  average  sum  by  this  number. 
"This  value  is  dependent  upon  the  gain  of" 

“the  device  and  will  probably  need  adjustment" 
“during  run.  Start  with  a  value  of  about  4003" 


PRINT  Pt 


ItVl1 


'i!i 


i 

s 

3413 

L v  1 :  INPUT  Ilvi 

« 

3423 

1 

V 

%• 

3433 

I 

yj 

3443 

!  set  up  the  plot 

1 

3453 

! 

1 

34S3 

! 

5 

T  4  -»A4 

j-mj 

Again : CALL  Plot_setup( Np lot ,S i tel .Maxrng , 3cale_y .Exrng , 2one$ ) 

& 

34S3 

Nrec"3 

5 

3433 

FOR  1-1  TO  233 

k?a 

3533 

ct; I >-3 

1 

3513 

NEXT  I 

1 

7T1/4 

JJLU 

Y-3 

3533 

Npoint 1 *0 

3543 

OUTPUT  KBD; *L“  1 

x 

3553 

1 

fl 

3553 

! 

i 

3573 

Sync:  !  synchronize  data  collection  with  clock 

3533 

T 1 -INT ( TIMEQATE  MOD  86433) 

i 

3533 

IF  T 1 <T3  THEN  T0-T3-8S400 

3533 

IF  T1-T0C1  THEN  GOTO  Sync 

* 

35 1  3 

T0-T1 

3523 

! 

g 

3633 

1 

1 

3543 

1  Data  collection  and  reduction 

5 

3653 

1 

•r 

V 

3553 

! 

<w 

3573 

Read_sig:  !  read  the  A-0  converter 

s 

3683 

CALL  Read_ad200(Dat<  » ) .Counts .TimeS ,0f f_ave  ) 

i 

3533 

Npoint"INT<  T1  MOD  3600  MOD  Nplot) 

! 

3733 

Npoint 2*Npoint -Npoint 1 

ti 

3713 

Nunl »Num-l 

3723 

Nun2=Cnt-Nun 

w* 

Jl  Jl] 

Ovd=*  1  .  /  ( Nun  ) 

0  t  43 

K-0 

1 

3753 

! 

1 

3763 

! 

V 

3773 

FOR  I-t  TO  Nun2  STEP  Nun 

X 

3783 

R«I*Rdist 

y 

3733 

IF  I  <Far  THEN  R-1 

7^77 
JO  O  O 

!<=K  +  1 

M 

7  i  /7 
j  w  •  O 

v  a=3 

I 

7  ^  ~  '"X 
w  *•  W 

FOR  J=I  TO  I ! Nunl 

7  “1  *7  ■> 

J  w  — 

X=Cat ( J  /-Offset 

3343 

Run_ave*Run_ave-t-X  »X 

s 

1853 

NEXT  J 

.■s 

A 

3863 

Ns=SQR( Run_ave»Qvd ) 

u 

■ 

CO 

-J 

Ci 

Point  <  K  )»R« < Ns-Corr ) 

S 

3833 

IF  K 1 > 1 0  THEN 

r* 

3893 

Ct(K  )-Ct(l<  H<  (Point (K  ) “2  ) «EXP(  2 »R»At ten  )  ) 

K 

3233 

Ct(K )-<Ct(K )»( (3.339E-8 )“2 ) »Er )/Zinp 

f- 

381  3 

ENO  IF 

3080 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4000 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4100 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4200 

4270 

4200 

4290 


IF  Point(K  )>ABS( 327S7  )  THEN  Point!  K >-32707 
D2(K+10)-INT(Point(K ) ) 

NEXT  I 

IF  K 1  > 1 0  THEN 
Y-Y+l 
END  IF 


Noise  correction  routine 


IF  Ns<<  Limit+S )  THEN 
Noi se*Noi se+Ns 
Kl-Kl+1 
IF  K 1  > I  0  THEN 

Corr*Noise/Kl 

Limit=*Corr 

ELSE 

Corr=*0 
END  IF 
END  IF 


Plotting  of  the  data 


Kmax»K 

REDIM  D2(Kmax+lS) 

! 

!  positioning  the  data  on  the  plot  by  tine  of  trac 
! 

T ime-< Npoint2/Nplot )»420 

Timedist*2<-(( Npoint2+Npoint  1  >/Nplot  )«420 

IF  Npoint1”0  THEN  Time*5 

IF  Npoint<E  THEN  Tinedist=G 

End_pt»Scale_y »2G0 

D2C  2  )  *  K  m  a  a 

D2<  3  >-INT( Temp  ) 

D2(  4  )*INT(  T me  ) 

D2( 5  >*INT( Timedist « 1 0  ) 

WINDOW  0,420,0,200 
GRAPHICS  ON 


!  set  the  vertical  increment 

! 

Inc»(  ( Num*Tminc*Rspd*End_pt  //Maxrng  ) 
02(0  )-INT< Inc  *  1 000 ) 


•TIN 


I 

I 


$ 


4400 

4410 

4420 

4430 

4440 

44S0 

44S0 

4470 

4480 

4490 

4900 

4510 

4520 

4530 

4540 

45SC 

4560 

4570 

4530 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

47C0 

4770 


I  compute  the  Intensity  of  the  return,  move  to  the  proper 
I  coordinates  and  plot  the  appropriate  colored  block 
! 

FOR  K-1  TO  Kmax 
Z-Point(K  >/Ilvl 
IF  ZM  THEN  Z-1 
IF  Z<0  THEN  Z-0 
AREA  INTENSITY  1,1,1 
Qis-INT(K»Inc  )  +  1 
MOUE  Timedist ,0iS 
RECTANGLE  Time  , Inc  .FILL 
NEXT  K 
! 

!  keep  an  account  of  the  trace  numbers  taken  on  the  plot 
I 

Nrec-Nrec+ I 
Rec_num-INT(  Nrec 1 
D2( 1 )«Rec_num 
Npoint I -Npoint 


save  routine 


if  function  key  is  set  then 
the  plot  will  be  saved 


IF  Cntrl-1  THEN 

ASSIGN  SFilel  TO  FilelS 
OUTPUT  BFilel ,Rec_num;02C  * ) 
END  IF 


graphics  dump  of  plot  after  15  minute  intervals 


IF  T imedist>415  THEN 
Prt_dmp :  PRINTER  IS  701 

PRINT  " 

PRINT  " 

PRINT  ” 

Plot  num=*P  lot  num  i  1 
DUMP  GRAPHICS  4701 
FOR  1-1  TO  Kmax 

Ct< I)-Ct< I )»( (Tempi273  5A2)*( t . / .0039  ! *K3»( 1  ./Y  ) 
CtC I )»Ct < I  )»< 1 ./(A*G  ) ) 

CtC I )-Ct< I  >»< 1 ./Pi  ) 

Ct( I )-Ct( I )/Pls_lng 
IF  Ct(I XI .0E-T00  THEN  Ct ( I  >=■ 1 . 0E- 1 00 
NEXT  I 


!  Ci  2  Computations  and  plots 


$ 


GCLEAR 

VIEWPORT  15,120,15,80 

WINDOW  0  ,300  ,0 .Exrng 

AXES  30  ,Exrng/ 1 5 ,0  ,0  ,30  ,Exrng/3 

CLIP  OFF 

CSIZE  2..S 

LORG  G 

FOR  M-0  TO  300  STEP  30 
MOVE  M  ,-Exrng/ 45 
LABEL  “1 . 0E“ ; < (M/30  )-t  0  > 

NEXT  M 

MOVE  140,-15 
CSIZE  4 


m 


LORG  8 

FOR  M-0  TO  Exrng  STEP  INT(Exrng/3) 
MOVE  O  ,M 
LABEL  M 
NEXT  M 
LOIR  PI/2 
LORG  G 

MOVE  -40  ,Exrng/2 
CSIZE  4 

LABEL  “RANGE  IN  METERS" 

LDIR  0 
LORG  4 

MOVE  1 50  ,Exrng+3 

LABEL  "TIME  AVERAGED  CT"2  “:Site$ 

CLIP  ON 

MOVE  300,0 

FOR  1-1  TO  Kmax 

Yp l-< Maxrng/Kmax  )*I 
Pnt-<LOG(Ct( I  )  )/2. 3025851  >  +  10 
IF  Pnt > 1 0  THEN  Pnt-!0 
IF  Pnt  <0  THEN  Pnt-0 
Pnt-Pnt  *30 
DRAW  Pnt.Ypl 
NEXT  I 


PRINT  “ 


PRINTER  IS  CRT 

GCLEAR 

Kl-10 

Lini t-Corr 
Noi se- 1 0  *Corr 


PM 


5380 

5393 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5473 

5480 

5490 

5500 

5510 

5523 

5530 

5540 

5550 

5553 

5570 

5580 

5590 

5503 

5610 

5623 

5533 

5540 

5650 

5560 

5670 

5630 

5630 

5700 

5710 

5720 

5733 

5740 

5730 


IF  CntrI-1  THEN 

ASSIGN  @Fi la  I  TO  * 
END  IF 
Cntrl=0 


if  the  plot  is  to  be  saved,  the  file  is  created 
and  named 


IF  Save_plt»l  THEN 

CALL  File_init<0i3c_address$  ,Nrec  ,Fi  le!  $  ) 
Save_pl t=0 
Cntr 1«1 
END  IF 


start  the  next  15  minute  plot 


GOTO  Again 
END  IF 
! 

I 

!  next  trace 
I 


GOTO  Sync 


toggle  set  if  plot  is  to  be  saved 


Svpt:  Save_plt»1 

GOTO  Read_sig 


completion  routine 


3730 

53C3 

5310 

6000 

6010 

6023 


END 

i 


SUBROUTINE  SECTION 


53 


tm 


0330  SUD  Freq_init ( FreqS ,N$  ) 

5043  ! 

5050  ! 

6050  f  SETUP  OF  THE  HP33I4A  FUNCTION  GENERATOR 
5073  ! 

5080  ! 

5090  '  Ap$*“AP" 

6100  FrqS-"FR* 

5110  Nm$»“NM“ 

5120  En$»'EN" 

5130  UoS*“U0" 

5140  HzS-'KZ* 

6150  INPUT  “FREQUENCY  DESIRED  (kHZ)  <5  RECOMMENDED  )•  .FreqS 

5150  INPUT  “AMPLITUDE  DESIRED  <U  .  . 1.SU  MAX)  *  ,Amp$ 

6170  IF  UAL( AmpS )> 1 . 5  THEN 

5180  Amp$““ 1.5“ 

5190  PRINT  "AMPLITUDE  OF  FUNCTION  GENERATOR  IS  1.5  U“ 

6200  ENO  IF 

6210  INPUT  “NUMBER  OF  CYCLES  PER  BURST  (INTEGER)  (100  RECOMMENDED >“  ,NS 

5220  OUTPUT  707i"MO3* 

6230  OUTPUT  707 ; “SR2“ 

5240  OUTPUT  707;Ap$&Amp$8,Uo$&Frq$8rFreq$&HzS&Nm$&N$SEn$ 

6250  SUBENO 
5250  I 
6270  ! 

5280  I 

6500  SUB  Init_ad200 
6510  1 

5520  I 

5530  !  INITIALIZATION  OF  THE  A-D  CONVERTER 

5540  ! 

5550  ! 

65S0  Ad_sel_code“! 7 

5573  Oummy*REAOIQ<  Ad_5al_c ode ,3 ) 

5530  WRITEIQ  Ad_3el_ccde ,0 ;0 

6330  CONTROL  Ad_sei_code ,0; t 

5500  SUBEND 
5513  ! 

5520  1 

5530  ! 


54 


7000 

7010 

7020 

7030 

7040 

7050 

70G0 

7070 

7080 

7030 

7100 

71  10 

7120 

7130 

7140 

7150 

7150 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7250 

7270 

7280 

7283 

n 

i  \J 

7313 


SUB  Set_t ine(  Zone* > 

! 

! 

!  SET  THE  TIME  DATE  RECORDER 
I 

I 

PRINT  “WHAT  TIME  REFERENCE  ARE  YOU  USING?  INPUT:’ 

PRINT  *  1  FOR  UNIVERSAL  TIME' 

PRINT  ’  2  FOR  LOCAL  TIME’ 

PRINT  ’  3  FOR  YOUR  OWN  CLASSIFICATION’ 

INPUT  K 
IF  K-2  THEN 

Zone*-* (LOCAL)*  ' 

ELSE 

IF  K-3  THEN 

INPUT  ’WHAT  IS  YOUR  TIME  REFERENCE’ .Zone* 

ELSE 

ZoneS-’(UTC)’ 

END  IF 
END  IF 

IF  TIMEDATE<OATE( '14  AUG  1984’)  THEN 
INPUT  ’ENTER  ”DD  MMM  YYYY’  *  *  .Date® 

INPUT  ’ENTER  ”HR:MIN:SC”  *  .Time* 

SET  TIMEDATE  DATE* Date* >+TIME< Time* ) 

PRINT  OATE*< TIMEDATE  )  .TIMES* TIMEDATE > 
Tatart-TIMEDATE 
TO-Tatart  MOD  8S400 
ENO  IF 

SUBENO 


55 


vC>: 


7500 

7SI0 

7520 

7530 

7540 

7550 

7SG0 

7570 

7580 

7590 

7S00 

76 1 0 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

77G0 

7770 

7780 

7790 

7800 

7310 

i  o<-  d 

7333 

n  4  'y 
I  u  **  v/ 

7353 

7360 

7370 


SUB  Read_ad200C INTEGER  0at< *  )  BUFFER .Counts .Time* , Of f_ave  ) 


INFOTEK  A-0  ROUTINE  SET  UP  FOR  EXTERNAL  TRIGGER 


Ad_ael_code-l 7 
I 

! INITIALIZATION  OF  THE  A-0  CONVERTER 
! 

OUTPUT  Ad_3el_codos “RESET” 

OUTPUT  Ad_sel_code ; *  INTERNAL' , “COUNT  "iCount* , "HOLDON* 
OUTPUT  Ad_ael_code ; “DELAYON’* , ‘SELECT  tal  end*  , Tine* 
OUTPUT  Ad_ael_code; "STATUS" 

ENTER  Ad_sel_code ;Resp* 

IF  Reap*-* - *  THEN 

ASSIGN  0Ad200  TO  Ad_aei_code;UQRD 
! 

(  triggering  of  the  HP3341A 
! 

IF  Of f_ave-0  THEN 
TRIGGER  707 
END  IF 

ASSIGN  @Buf  TO  BUFFER  Dat(») 

TRANSFER  8Ad200  TO  0Buf;WAIT 
OUTPUT  Ad_ael_code;“  * 

OUTPUT  Ad_ael_code! "STATUS" 

ENTER  Ad_ael_code ;Reap* 

IF  ReapSO* - -  THEN 

PRINT  "ERROR-  "sReap* 

END  IF 

ELSE 

PRINT  “ERROR  DURING  INITIALIZATION  -“ iReap* 

ENO  IF 
3UDEND 

I 

! 

! 


56 


8800 

8010 

8020 

8030 

8040 

8050 

8060 

3070 

8080 

8090 

8100 

8110 

8120 

8130 

8140 

8150 

8160 

8170 

8180 

8150 

8200 

8210 

8220 


SUB  Plot_setup(  Nplot  .Site*  .Maxrng  ,Scale_y  .Exrng  , Zone®  ) 
I 
! 

!  SET-UP  OF  THE  TIME  PLOT  ON  THE  CRT 


Scale_y-Maxrng/Exrng 
GRAPHICS  ON 
LINE  TYPE  I 
VIEWPORT  15,120,15,80 
WINDOW  0  .Nplot  ,0  .Exrng 

AXES  Nplot /I  5 .Exrng/ 15 ,0  ,0  ,Nplot/3  ,Exrng/3 
CLIP  OFF 
CSIZE  4, .6 
LORG  6 

T1-TIMEOATE  MOO  86400 
Hrs-Tl  QIV  3600 
T2-T1  MOO  3600 
Min-T2  DIU  60 
Qtrhr-Min  OIV  IS 

FOR  M-0  TO  Nplot  STEP  TNT(Nplot/3) 

MOVE  M,-Exrng/45 
Qtrmin-Qtrhr * 15+<  P1»3/Nplot  )»5 


8230 

8240 

8250 

8260 

8270 

8230 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

83G0 

8370 


IF  Qtrmin-60  THEN 
Qtrmin-0 
Hrs-Hrs+1 
END  IF 

LABEL  USING  *00  ,A  ,ZZ* tHraj  * : * iQtrmin 
NEXT  M 

MOVE  Nplot/2,-15 
LABEL  ‘TIME  "SrZone* 


I  LABEL  ORDINATE 
! 


LORG  8 

FOR  M-0  TO  Exrng  STEP  INT(Exrng/3> 
MOVE  0  ,M 
LABEL  M 


0300 

NEXT 

M 

0300 

LOIR 

PI/2 

{J-tiJU 

LORG 

G 

8410 

MOVE 

-Np lot/7  ,Exrng/2 

8420 

LABEL  “RANGE  IN  METER 

8430  ! 

8440  ! 

TITLE 

8450  ! 

8460 

LOIR 

0 

8470 

LORG 

4 

8400 

MOVE 

Nplot/2  ,Exrng  +  3 

{Sites’ 


8490 

8900 

8510 

8520 

8530 

8540 

9000 

9010 

9020 

9030 

9040 

90S0 

9050 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

91G0 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

9240 

9250 

9250 

9270 

9280 

9290 

9300 

9310 

9338 

9340 

9350 

9350 

9370 

3380 

9390 

9400 

9410 


LABEL  DATES  ( TIMEDATE 1  s  "■ 
CLIP  ON 
SUBEND 


SUB  Attenuation*  Freq  .Temp  .Atten  ) 


!  This  subprogram  calculates  the  attenuation  of  the 
!  sound  in  air  based  upon  equations  in  Neff  1975 
!  ( source  of  subroutines  Thesis  of  R.  Fuller) 

I 

1  Variables 

I 

1  Atom_pres  -  input  atmospheric  pressure  in  mb 

! 

!  Atten  -  attenuation  coefficient  of  acoustic  wave 
! 

!  Att_max  -  Variable  in  program.  It  ia  the  attenuation 
!  at  the  frequency  of  the  maximum  attenuation 

I  for  the  given  input  conditions. 

1 

!  F  -  Ratio  of  the  frequency  to  frequency  at  maximum 
I  attenuation. 

! 

!  Fmax  -  Frequency  of  the  maximum  attenuation 
I 

1  H  -  variable  used  in  the  integration  of  excess  attenuation 

! 

!  Pstar  -  variable  used  in  intermediate  calculations 
! 

!  Tstar  -  variable  used  as  an  intermediate  in  calculation 
I  of  attenuation. 

! 

!  Uat_pres  -  Atmospheric  water  pressure  in  mb. 

! 

INPUT  "Enter  the  atmospheric  pressure  in  mb" ,Atom_pres 
INPUT  "Enter  the  atmospheric  water  pressure  in  mb" ,Wat_pres 
H- 1 00»Wat_pres/Atom_pres 
Tstar-* 1 .  3*Temp  +  49Z )/51 3 
Pstar-Atom  _pres/1014 

Fmax-< 1 0+SS00*H+44400*H»H ) *P star /Tstar * .  3 

At t_max« .0078 *Fmax* Tstar " * -2.5 >*EXP( 7. 77*  * 1  - 1 /Tstar ) 5 

F-Freq/Fmax 

Atten-* At  t_max/304 .  3  )*<  <\  f8*F  rZ+C2*F»F/<  l+F*F>  )"2  )".S 

Atten-< Atten+1 .74E-10*Freq*Freq  )/4.35 

SUBENO 
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9430  I 

9500  SUB  File_init ( Oisc_addressS ,Nrac ,Fi le 1  $ ) 

9510  1 

9520  ! 

9530  I  CREATE  THE  STORAGE  FILE  ON  THE  DISC  FOR 

9540  !  THE  REDUCED  DATA 

9550  I 
9560  I 

9570  INPUT  “ENTER  THE  REDUCED  OATA  OUTPUT  FILENAME  “  File* 
9580  Fi  le!S«FileS8r0i3C_addras5* 

9S90  CREATE  BOAT  Filel S , 180 ,400  • 

9600  ASSIGN  SFilel  TO  Filel* 

9610  ! 

9620  SUBENO 


SPEAKER  AND  ARRAY  ANALYSIS 


In  the  design  of  our  echosounder,  it  was  determined 
that  a  rapid  decay  time  was  required  to  obtain  accurate 
short  range  information.  Weight  restrictions  involved  with 
equipment  transportation  plus  maximum  response  at  high 
frequencies  led  to  our  decision  to  use  piezo  ceramic 
speakers.  The  Motorola  KSN  1005A  speaker  was  selected 
based  on  these  requirements  and  the  specifications  charted 
in  the  Motorola  catalog  [Ref.  16]  and  reproduced  below  in 
Figures  23  through  25. 


Nominal  Power,  Impedance,  and  Distortion  Ratings 


FREQUENCY  (Hz) 

300  500  1000  2000  5000  10000  20000  40000 

Fig.  23.  Speaker  Ratings 


The  frequency  response  charts  indicate  that  a  maximum 
response  for  our  speakers  occurs  at  a  resonant  frequency  of 
5000  Hertz.  This  frequency  was  used  as  the  baseline  from 
which  all  our  measurements  are  made. 

Using  the  speakers  in  an  anechoic  chamber,  the  average 
e- 1  voltage  decay  time  was  measured  to  be  approximately  900 
usee  (Fig.  26).  This  decay  time  translates  into  a  sound 
propagation  distance  of  just  over  15.0  centimeters  (at  STP) 
from  the  speakers.  Considering  our  requirements,  this 
speaker  is  ideally  suited  to  serve  our  purpose. 


B.  ACOUSTIC  ARRAY 


Our  next  consideration  was  the  echosounder  array 
pattern.  Ideally,  the  acoustic  sources  should  be  placed 
exactly  one  half  wavelength  apart.  At  a  frequency  of  5000 
Hertz,  this  would  require  spacings  of  3.4  centimeters  (at 
STP)  which  is  physically  impossible  for  the  speakers  we 
have  chosen.  The  closest  possible  spacing  is  7.62 
centimeters  between  sources  after  shaving  off  the  flange  of 
the  horn  (Fig.  25). 

From  Kinsler,  et  al .  [Ref.  14],  the  equation  for  the 
directionality  factor  of  a  simple  line  array  is  derived  as: 

N 

sin  (  -  kd  sin  8  ) 

2 


1 

sin  (  -  kd  sin  8  ) 

2 

where 

k  is  the  wave  number  ( 2x/ a  ), 
d  is  the  distance  between  sources, 

N  is  the  number  of  sources,  and 

8  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest. 

However,  this  equation  assumes  simple  point  sources  which 
does  not  adequately  describe  our  speakers.  It  was  necessary 
to  couple  this  equation  to  the  directional  factor  for  a 


1 

H  ( N ,  8  )  = 

N 


s 

$ 


piston  source  which  is  also  identified  in  Kinsler,  et  al . 
[Ref.  14:  p.  108]  as: 


D(0) 


2  Ji  (  ka  sin  (  8  )  ) 
ka  sin  (  9  ) 


where 

k  is  still  the  wave  number, 

a  is  the  radius  of  the  piston  source, 

9  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest,  and 

Ji  (  x  )  is  a  Bessel  function  of  Order  1  with  an 
argument  x . 

These  two  equations  were  combined  together  to  form  the 
equation  of  directionality  for  a  linear  array  of  piston 
sources,  L(N,8),  by  simple  multiplication. 

L(N,9)  =  H(N, 0 )  •  D(0) 

This  equation  was  incorporated  into  a  computer  program, 
which  may  be  found  in  Reference  3,  and  used  to  generate 
linear  array  patterns  for  various  numbers  of  piston  sources 
(Figs.  27  through  31). 

It  was  concluded  from  these  plots  that  five  linear 
elements  would  provide  the  best  combination  of  forward 
directionality,  sidelobe  suppression,  physical  size  and  a 
relatively  low  cost.  Then,  in  order  to  enhance  both  array 
efficiency  and  symmetry,  we  settled  on  a  five  by  five 


element  array  design  with  7.26  centimeter  spacing  between 


decibels 


speakers  in  both  the  vertical  and  horizontal  directions. 


After  verifying  the  manufacturers  polarity  designation 
for  35  speakers,  we  obtained  Lissajous  plots  for  each 
individual  speaker  (Fig.  32).  Based  on  the  speaker's  output 
to  input  voltage  ratios  as  illustrated  by  these  plots,  we 
were  able  to  rank  all  our  speakers  by  signal  efficiency. 
It  was  based  upon  this  criteria  that  we  selected  the  25 
most  efficient  speakers  for  the  array,  placing  the  best 
speakers  at  the  center  and  subsequent  ranking  speakers 
further  toward  the  sides  and  corners. 


Fig.  32.  Lissajous  Plot 


The  25  selected  speakers  were  r'.ounted  in  a  five  by  five 
planar  array  on  a  balsa  wood  insulated  bilayered  sheet 
metal  board.  After  wiring  all  the  speakers  in  parallel,  we 
surrounded  all  the  electrical  connections  and  speaker  backs 
with  two  3.0  centimeter  layers  of  foam  insulation 
sandwiching  a  1.0  millimeter  lead  sheet.  Then  the  entire 
array  mounting  was  enclosed  in  a  44  by  44  by  5  centimeter 
sheet  metal  box.  This  design  was  chosen  to  suppress 
virtually  all  acoustic  energy  propagating  out  the  rear 
hemisphere  of  the  array,  while  shielding  the  array  from  any 
external  electrical  interference  (Figs.  33  and  34). 


Fig.  33.  Array  Photo 
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ENCLOSURE  DESIGN 


: 


I 


Acoustic  echosounding  has  proven  to  be  an  extremely 
useful  technique  for  probing  and  analyzing  the  lower 
atmosphere.  In  order  to  most  efficiently  utilize  the 
acoustic  waves  transmitted  and  later  received  by  this 
remote  sensing  method,  it  is  essential  to  have  an  efficient 
antenna  with  highly  directive  beams  and  strongly  suppressed 
sidelobes.  Antenna  design  becomes  increasingly  more 
important  in  a  noisy  environment  where  noise  pollution 
within  the  sidelobes  may  dominate  the  desired  signal  within 
the  main  lobe.  Hall  and  Wescott  [Ref.  12]  showed  that 
sidelobe  suppression  improved  with  higher  frequencies. 
Their  studies  showed  that  the  measured  90  degree  sidelobe 
suppression  ranged  from  38  dB  at  1  KHz  to  50  dB  at  5  KHz. 
Furthermore,  any  significant  improvement  in  sidelobe 


ft 


suppression  could  only  be  obtained  by  surrounding  the 
antenna  with  an  acoustic  energy  absorbing  cuff  or  shroud. 

In  an  effort  to  maximize  our  antenna  main  lobe  to 
sidelobe  power  ratio,  we  intend  to  operate  only  at  high 
frequencies  as  discussed  in  the  previous  section. 
Additionally ,  we  have  designed  an  acoustic  energy  absorbing 
enclosure.  Many  designs  were  considered  based  upon 
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previous  research  in  the  field  of  echosounding  [Refs.  17 
through  20],  In  addition,  we  obtained  the  actual  acoustic 
beam  patterns  for  our  array  using  a  computer  program 
written  by  LCDR  Butler  [Ref.  21]  which  we  modified  for  our 
purposes.  This  modified  version  of  LCDR  Butler’s  program 
may  be  found  in  Reference  3.  By  rotating  the  array  in  an 
anechoic  chamber,  we  were  able  to  produce  highly  accurate 
polar  plots  of  the  array  beam  patterns  (Figs.  35  and  36). 


5X5  ACOUSTIC  ARRAY  BEAM  PATTERN 


R-3.60  M  MIKE  VT-  8.2465  INPUT-5.0  V  FREQ-5. 0  KHZ 
Fig.  35.  Polar  Plot  of  5  X  5  Array  Beam  Pattern 


71 


These  polar  plots  and  the  actual  coordinates 
corresponding  to  the  individual  data  points  confirmed  the 
computer  prediction  for  the  five  element  line  array  (Fig. 
29).  Additional  polar  plots  obtained  by  varying  both  the 
input  array  voltage  and  range  between  array  and  microphone 
further  support  our  claim  that  for  a  5  KHz  carrier 
frequency,  the  main  lobe  is  confined  to  a  divergence  angle 
of  20  degrees.  Since  it  is  our  aim  to  suppress  all 
sidelobes  and  utilize  solely  the  main  lobe,  we  chose  not  to 
taper  our  enclosure  as  most  previous  researchers  had. 
Rather  we  designed  the  enclosure  based  upon  the  dimensions 
of  the  array  itself  and  the  acoustic  beams  it  generated. 

Plywood  was  used  for  the  construction  of  the  enclosure 
and  provided  not  only  a  rigid,  inexpensive  framework,  but 
also  proved  to  greatly  attenuate  external  noise 
interference.  Anticipating  all  kinds  of  weather  conditions 
during  data  collection,  the  plywood  enclosure  was  first 
waterproofed  with  four  coats  of  marine  varnish.  Grooved 
joints,  caulking  and  weather  stripping  were  also  design 
considerations . 

A  millimeter  layer  of  lead  can  suppress  an  acoustic 
signal  as  much  as  40  dB  (Figs.  37  and  38).  About  7.0 


centimeters  of  corrugated,  egg-carton  design  foam  can 
suppress  a  signal  another  3  to  4  dB  (Figs.  37  and  39). 
Together  they  make  an  extremely  efficient  absorbing 


3ThRT :  4  O00  Hz 

BH:  1?.09~  Hi 

STOP:  i  000  Hi 

?000  Hi 

Yi  -!■?.  ^3  ctBS' 

Fig.  37.  Response  Reference,  No  Insulation 


RANGE!  -21  dBV 


STATUS:  PAUSED 
RMS: 23 


-101  i _ ; _ _ _ « _ i _ - - 

START:  4  000  M2  BU:  19.097  Hz  STOP:  6  000  Hz 

Xl  5000  Hz  _ Yi  -22. 78  JBV _ __ 


Fig.  39.  Signal  Suppression  by  Foam 

Lead-lined  absorbing  foam  was  glued  to  all  inner 
surfaces  of  the  enclosure  with  two  1.0  millimeter  layers  of 
lead  overlapping  at  all  corners.  Strong  aluminum  brackets 
were  used  to  connect  the  four  side  panels  to  each  other  as 
well  as  to  the  enclosure  base  (Figs.  40  through  42). 
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Fig.  40.  Enclosure  Photo,  Fully  Assembled 


APPENDIX  D 

ECHOSOUNDER  OUTPUT 

Sixteen  echosounder  output  traces  were  included  to 
exhibit  typical  atmospheric  activity.  Many  of  the  plots, 
such  as  the  14:45,  15:30,  and  the  17:30  of  26  April  and  the 
10:00,  11:15,  13:00,  14:15,  15:15,  16:15,  and  the  17:15  of 
27  April,  have  convective  plumes  which  are  prevalent 
whenever  a  heat  flux  between  the  surface  and  atmosphere 
exists.  The  plots  of  18:15  and  19:15  on  the  26th  of  April 
and  of  18:00,  18:45  and  20:00  on  the  27th  of  April  clearly 
show  the  passage  of  the  neutral  event  which  is  encountered 
when  the  atmospheric  and  surface  temperature  difference 
becomes  negligible.  Finally,  the  plot  of  16:30  on  26  April 
can  be  associated  with  strong  winds  which  exhibit  a 
somewhat  uniform  return  for  all  altitudes  across  the  entire 
15  minute  interval  of  the  trace. 

The  horizontal  axis  labels  of  the  Ct 2  plots  are  small 
and  may  be  hard  to  read.  This  axis  is  a  log  scale 
beginning  with  1»10-10  at  the  left  side  of  the  plot  and 
ending  with  1*10°  at  the  right  side  of  the  plot.  Each  tick 
mark  moving  left  to  right  along  this  axis  represents  an 
integer  increase  in  the  exponent  of  10. 
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